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Abstract: We analyze the recently published Fermi-LAT diffuse gamma-ray measurements in the 
context of leptonically annihilating or decaying dark matter (DM) with the aim to explain simul- 
taneously the isotropic diffuse gamma-ray and the PAMELA, Fermi and HESS (PFH) anomalous 
data. Five different DM annihilation/decay channels 2e, 2//, 2t, 4e, or 4/i (the latter two via 
an intermediate light particle </>) are generated with PYTHIA. We calculate both the Galactic and 
extragalactic prompt and inverse Compton (IC) contributions to the resulting gamma-ray spectra. 
To find the Galactic IC spectra we use the interstellar radiation field model from the latest release of 
GALPROP. For the extragalactic signal we show that the amplitude of the prompt gamma-emission 
is very sensitive to the assumed model for the extragalactic background light. For our Galaxy we 
use the Einasto, NEW and cored isothermal DM density profiles and include the effects of DM sub- 
structure assuming a simple subhalo model. Our calculations show that for the annihilating DM 
the extragalactic gamma-ray signal can dominate only if rather extreme power-law concentration- 
mass relation C{M) is used, while more realistic C{M) relations make the extragalactic component 
comparable or subdominant to the Galactic signal. For the decaying DM the Galactic signal always 
exceeds the extragalactic one. In the case of annihilating DM the PFH favored parameters can be 
ruled out by gamma- ray constraints only if power-law C{M) relation is assumed. For DM decaying 
into 2/z or 4/i the PFH favored DM parameters are not in conflict with the gamma-ray data. We 
find that, due to the (almost) featureless Galactic IC spectrum and the DM halo substructure, 
annihilating DM may give a good simultaneous fit to the isotropic diffuse gamma-ray and to the 
PFH data without being in clear conflict with the other Fermi-LAT gamma-ray measurements. 
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1. Introduction 

During the last few years several experiments have shown an anomalous excesses in the 
cosmic electron and positron spectra. The PAMELA satellite has observed a steep rise of 
positron fraction +6"*") at energies above 10 GeV with no significant excess in the 

cosmic antiproton flux [1, 2]. The Fermi satellite and the HESS atmospheric Cherenkov 
telescope have measured an excess of high-energy (e~ + e"*") flux with a cut-off of around 
800 GeV [3, 4]. The ATIC and PPB-BETS balloon measurements indicate a similar excess 
[5, 6]. Most excitingly, the excess might originate from the annihilation or decay of the 
dark matter (DM) particles. The nature of those signatures requires the properties of 
DM to deviate strongly from the standard freeze-out predictions. The thermally averaged 
DM annihilation cross-section {ax^) has to be boosted some orders of magnitude over the 
standard freeze-out value (crAi^)std — 3 x 10~^^ cm^s~^, which might be achieved, e.g., 
through the Sommerfeld effect [7] (see, e.g., [8, 9, 10] for the related phenomenological 
studies) or through the Breit-Wigner resonant enhancement [11, 12, 13, 14]. On the other 
hand, the decaying DM [15, 16, 17, 18, 19] can explain the excess independent of the freeze- 
out constraints. In both cases, the annihilation or decay of DM should favorably occur only 
through the leptonic channels [20, 8, 21, 22, 23], as no excess in the hadronic channels has 
been observed. Alternatively, the excess of can potentially be explained by modifying 
or adding astrophysical sources, e.g. pulsars [24, 25, 26, 27]. 

The high energy leptons of the DM annihilation/decay are inevitably accompanied 
by the gamma-rays due to the final state radiation of charged leptons and decays of sub- 
products ("prompt gamma-rays") and due to the upscattered background photons from 
the inverse Compton (IC) scattering ("IC gamma-rays"). Thus, the observed gamma-ray 
fluxes strongly constrain the above mentioned cosmic ray anomalies from DM annihila- 
tion/decay. The strongest gamma-ray constraints should arise from the observations of the 
Galactic center (GC) [28, 29], as the density of DM is very high and it is relatively close 
to us. Those analyses take into account both the prompt and IC gamma-ray contribu- 
tions [30, 31, 32, 33, 34, 35]. On the other hand, the GC is densely populated by different 
astrophysical objects, which contaminate the gamma-ray signal and introduce significant 
uncertainty in the derived constraints. The Galactic gamma-ray signal of DM annihila- 
tion/decay at higher latitudes is considerably weaker than the signal from the GC. Despite 
being considerably weaker, the suppressed contamination by the Galactic astrophysical 
sources partially compensates the weakness of the signal. In addition, adding the Galactic 
DM substructure into the picture [36, 37, 38, 39, 40] may significantly change both the 
magnitude and the morphology of the induced gamma-ray signal as the diffuse signal from 
DM subhalos can be essentially isotropic. The Galactic gamma-ray signal of DM annihi- 
lation/decay at higher latitudes can also include a considerable extragalactic contribution. 
It can originate from different sources: DM annihilation/decay in cosmological distances, 
active galactic nuclei (AGN), structure formation shocks, starburst galaxies, etc. 

Fermi-LAT collaboration has also derived constraints on the annihilating DM proper- 
ties from the Galactic dwarf spheroidal galaxies [41] (see [42, 43] and references therein for 
the other recent studies) . Non-detection of gamma-ray signal towards several nearby galaxy 
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clusters has a potential to severely constrain the annihilating DM models [44, 45] . However, 
the present Fermi results [46] still allow for the DM parameters giving the best fit to the 
PAMELA, Fermi, and HESS (PFH) anomalies even if the smallest DM substructures 
are present. In addition to gamma-rays, the energetic leptons from DM decay/annihilation 
produce other accompanying signatures: synchrotron radiation [30, 47, 48], and neutrinos 
from fj, OT T decay [49, 50, 51, 52, 53]. The gamma-rays ionize and heat the intergalactic 
gas. The additional electrons released in the ionization process change the scattering op- 
tical depth of the CMB photons [54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65]. Also, the 
energetic released in the annihilation/decay process induce the nonthermal component 
to the Sunyaev-Zel'dovich effect in galaxy clusters [66, 67, 68]. The energy injection from 
the DM sector modifies the accurately calculable^ standard cosmic recombination process 
(e.g. [71, 54, 72]) and the Big Bang Nucleosynthesis (e.g. [73, 74, 75]). 

As in, e.g. Regis &: Ullio [76], in this study we focus on the gamma-ray signal of DM an- 
nihilation/decay at higher Galactic latitudes. Fermi-LAT collaboration published recently 
their measurement of the diffuse gamma-ray emission at several Galactic regions and their 
estimation of the "extragalactic diffuse gamma-ray background" [77]. The published spec- 
tra are more constraining and have considerably smaller error bars than the older analogous 
spectra from the EGRET experiment [78, 79]. The main aim of this paper is to find out 
whether the new Fermi-LAT diffuse gamma-ray measurements can rule out or, instead, 
to support the PAMELA, Fermi, and HESS favored models of DM annihilation/decay as 
given in [34]^. 

Our study is "model independent": we assume that DM particles annihilate/decay 
into Standard Model charged leptons, £~ , £~^; £ = e, ji, r, extending the analyses presented 
in series of works, e.g. [8, 9, 33, 34, 81, 82, 83, 84]. We analyze five channels presented 
in Table 1. Our general formalism for the extragalactic gamma-rays is presented in our 
previous work [60] . The new ingredient in this paper is to include the effect of the absorption 
of extragalactic gamma-rays due to the pair production on extragalactic background light'^. 
For that we test different models for the extragalactic background light and show that the 
signal of prompt gamma-rays from the DM annihilation/decay can be significantly reduced 
in some cases. For calculation of the Galactic IC spectra we use the interstellar radiation 
field (ISRF) model from the latest release of GALPROP [86]. We show that inclusion of the 
realistic ISRF including CMB, infrared radiation from dust, and stellar light backgrounds 
are complicated enough to smear out any spectral feature of the individual components 
and the resulting spectrum is essentially power-law-like. To model the Galactic DM halo 
we use three density profiles: Einasto, NEW and cored isothermal. We first perform our 
analyses assuming that our Galaxy consist of one structureless DM halo with the given 
density profile. This is probably unrealistic but often used approximation. After that we 
repeat our study using a more realistic model including DM subhalos. The technical details 
of our analyses are collected in Appendices A, B, C. 

^See, e.g., [69, 70] for the latest developments. 
^For an earlier analysis, similar to [34], see [80]. 

^See [85], where the authors have included the background light while discussing the prospects for 
detecting the decaying DM with the Fermi-LAT. 
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After Fermi collaboration has made their data publicly available, several papers have 
appeared that estimate the constraints on annihilating/decaying DM properties [82, 83, 
87, 88, 89, 90]. The first two papers [82, 83] use the preliminary Fermi-LAT data. In 
both of these studies the authors have neglected the potential extragalactic contribution 
while deriving the constraints for the annihilating DM. In [83] the approximate ISRF of 
our Galaxy is used. However, we show in this paper that the use of more realistic ISRF 
can significantly change the (line of sight) integrated IC spectra and in some cases the 
extragalactic diffuse gamma-ray signal may dominate over the Galactic one. Abazajian et 
al. [89] use the results of [77] but consider only the prompt gamma-ray spectra to constrain 
their DM models. This is clearly unsatisfactory because, as we show in agreement with 
[90], in general the bounds are dictated by the IC contribution to the diffuse gamma- 
ray spectrum while the prompt contribution plays just a subdominant role. The two-peak 
structure of the diffuse gamma-ray spectrum in DM annihilations into 2// channel presented 
in Fig. 4 of Abdo et al. [90] also shows that those authors do not use realistic ISRF to 
calculate the IC spectrum but approximate it with the CMB component. As we show in 
this paper, this approximation does not qualitatively change the bounds derived in that 
paper but may qualitatively affect the whole interpretation of the Fermi-LAT experimental 
results. In addition, both Abazajian et al. [89] and Abdo et al. [90] study only the 
annihilating DM. In this work, similarly to Refs. [43, 91], we study the difference between 
annihilating and decaying DM gamma-ray signals. 

The most important new result of this work is that the inclusion of realistic IRSF 
and Galactic DM subhalo structure into the analyses opens up a new interpretation of the 
Fermi-LAT isotropic diffuse gamma-ray measurement. We find that (i) the leptonically 
annihilating DM models can give good fits to the Fermi-LAT isotropic diffuse gamma- 
ray data if the extragalactic contribution is subdominant compared to the Galactic one; 
(ii) the best-fit regions of the preferred DM mass and annihilation cross section has an 
overlap with the best-fit regions of the PFH anomaly [34] if the boost factor from the 
Galactic DM substructure is i?sub ^ C(10); (iii) in the case of decaying DM the fits to the 
isotropic diffuse data are worse and there is no overlap with the best-fit regions of the PFH 
anomaly. This result implies that a significant fraction of the isotropic diffuse gamma-ray 
signal observed by Fermi-LAT may actually be of the Galactic origin. A generic issue that 
for a good fit the central regions of our Galaxy should give less prominent contribution to 
the DM annihilation signal than generally expected is (partially) solved with the DM halo 
substructure. The simultaneous fitting of PFH electron/positron and Fermi-LAT diffuse 
gamma ray data gives a preference to the DM models with Sommerfeld enhancement due 
to a new light intermediate particle (p [9] (see Table 1). If the intermediate particle (p is 
long-lived [92], a scenario not considered in this work, their long diffusion length further 
smears any DM annihilation signal and helps to resolve potential conflicts with the diffuse 
gamma-ray data measurements from the central Galactic regions. 

Our paper is organized as follows. In Section 2 we describe the energy and particle 
input from the DM annihilation/decay. In Section 3 we solve the radiative transfer equation 
for the Galaxy and the intergalactic medium. Section 4 presents the constraints as inferred 
from the Fermi-LAT diffuse gamma-ray data. Section 5 discusses the possibility of fitting 
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the models to the Fermi-LAT isotropic diffuse data. Our summary and final discussion is 
given in Section 6. 

2. Energy input from DM annihilation or decay 

In this paper we consider the DM annihilation and decay to Standard Model charged 
leptons presented in Table 1. In addition to the three two-lepton channels, which were also 
investigated in our previous work [60], we have now included two additional four-lepton 
annihilation and decay channels. In our calculations the intermediate particle (j) has a 
mass 7710 = 1 GeV, which forbids (j) to decay into r leptons, and thus we do not consider 
a 4t channel. The decay of (p (assumed to be prompt) and the final state distributions 
of e^, photons and u from the primary two-body states produced in DM annihilations 
are computed using PYTHIA 8.1 Monte Carlo tool [93]. Our DM annihilation scenarios 
for 2e, 2/i and 2r are the same as presented in [8]. The energy partition between three 
main products of the annihilation/decay process: (i) electrons/positrons, (ii) photons, and 
(iii) neutrinos/antineutrinos is given in Table 2. In channels involving unstable muons 
and taus a large fraction of energy is carried away by neutrinos which is lost for the 
creation processes of the gamma-ray background. The energetic electrons and positrons 
instantaneously interact with the ambient photon fields and create a soft part of the output 
gamma-ray spectrum via the inverse Compton (IC) mechanism. The most energetic part 
of the gamma spectrum consists of the prompt (or direct) photons originating directly from 
the annihilation/decay process. Table 2 shows that the 2r channel has a particularly high 
fraction of energetic prompt photons. Finally, it is important to mention that the spectra 
from annihilation and decay are related to each other in a simple manner: annihilation of 
two DM particles with a mass can be seen as a decay of a particle with mass 2m£,M, 

i.e. all of the decay spectra are indistinguishable from the annihilation spectra for two 
times smaller particle mass. 



Annihilation 


Decay 


Notation 


2DM ^ e+ e- 


DM -^e+ + e- 


2e 


2DM fi+ + fi- 


DM fj,+ + fi- 


2 11 


2DM -^T+ + T- 


DM ^ T+ + r- 


2t 


2DM ^ 2(j), ^ e+ + e- 


DM 2(t), (j) ^ e+ + e- 


4e 


2DM 2^, (j) ^ fi+ + /i" 


DM ^ 24>, 4> ^ fi+ + i^r 


All 



Table 1: The annihilation and decay channels of DM and their notation in the paper. 



3. Gamma-ray spectra from DM annihilation/decay 

In this section we briefly describe the machinery used to calculate the gamma-ray spec- 
tra. We start by introducing the approach used for calculating the contribution from our 
Galaxy's DM halo, and in the second half of the section we focus on extragalactic part. 
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channel 


e± 


photons 


neutrinos 


2e 


~ 96 - 97 % 


~ 3 - 4 % 


0% 


2fi 


~ 34% 


-2-3% 


~ 63 - 64 % 


2t 


~ 16% 


~ 16-17% 


~ 67 - 68 % 


4e 


~ 98% 


-2% 


0% 


4/u 


~ 35% 


-1% 


~ 64% 



Table 2: Energy partition between e*, photons and neutrinos/antineutrinos for all the studied 
annihilation and decay channels assuming rriDM = 100 GeV - 10 TeV. 

3.1 Galactic halo contribution 

As described above, the energetic created in the annihilation/decay interact with the 
ambient photon fields through the IC mechanism, where the soft input photons are boosted 
to gamma-ray energy range. Inside the Galaxy the interstellar radiation field (ISRF) 
relevant for the IC mechanism has three main components: (i) stellar light, (ii) infrared 
light, which is the reprocessed form of the stellar radiation by the surrounding dust, (iii) 
cosmic microwave background (CMB) radiation. In our calculations we take the ISRF 
from the latest release of GALPROP^ [86]. Since our input can be very energetic (we 
consider mr>M up to 10 TeV) the IC calculations involving stellar and infrared components 
of the ISRF need to exploit the full energy dependent Klein-Nishina form for the cross- 
section. For our IC calculations we use the formalism in Blumenthal & Gould [94], and 
in Cirelli & Panel [33]^. We have neglected the effect of electron diffusion, i.e. our IC 
photons are created "on spot"^. Once the gamma photon is produced, whether directly 
from the annihilation/decay event or through the IC mechanism, it can propagate almost 
freely inside the Galaxy, i.e. inside the Galaxy at gamma-ray energy range we can safely 
neglect any absorption. Thus, to calculate the gamma-ray intensity towards some Galactic 
latitude b and longitude I one just has to sum up the emissivity Je [eV cm~'^s~^eV~^] along 
the line of sight 

lEib,l) = ^[ jEib,l,r)dr. (3.1) 

The ISRF from GALPROP code is given as a grid in Galactic cylindrical coordinates 
R and z, covering the ranges R = 0.25 — 20.25 kpc and z = —5 — 5 kpc with a total 
of 101 X 41 grid points. We calculate the emission coefficients Je for each of the grid 
points'' by looking first at annihilating DM models with a standard thermal production 
cross section (crA'^)std — 3 x 10~^^ cm^s~^ and fixing the DM density globally to Pq = 0.4 
GeV cm~^, which is a recently preferred value for the DM density in our local Galactic 
neighborhood [96, 97]. We perform this calculation for all five annihilation channels and for 

^http://galprop. stanford.edu/ 

^For the sake of completeness we have given a "quick recipe" for the IC calculations in Appendix A. 
Also, Ref. [66] presents the IC formalism. 

®For more rigorous treatment one should first solve the diffusion-loss equation for and then calculate 
the IC spectrum using this new electron distribution. 

''in reality we can discard half of the grid points due to the reflection symmetry of the ISRF with respect 
to the z = plane. 
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Figure 1: Galactic gamma-ray spectra (solid lines) for 2/z annihilation channel with Einasto density 
profile, niuM = 1 TeV and with annihilation cross-section (cr^u) = 60 x 3 x 10"^^ cm'^s"^ for 
10 < 6 < 20 (upper left panel), 20 < 6 < 60 (lower left panel), and 6 > 60 (upper right panel) 
regions. Spectrum of the Galactic isotropic component is given on the lower right-hand panel. The 
decomposition of the Galactic spectra into main components (in increasing order of energy): (i) 
IC from CMB, (ii) IC from infrared radiation, (iii) IC from stellar light, (iv) prompt emission, 
is also presented. Dashed and dotted lines represent extragalactic spectra for the Maccio et al. 
[95] C(M) relation with Mmin = 10"*^ Mq and for the power-law C{M) with Mmin = 10 ^Mq, 
respectively. The high-energy splitting of these curves corresponds to the three choices for the 
extragalactic ultraviolet (UV) background (from top to bottom): (i) "no UV", (ii) "realistic UV", 
(iii) "optimistic UV" . The points with errorbars give the Fermi-LAT measurements. 



20 logarithmically-spaced DM particle masses between 100 GeV and 10 TeV. Our emission 
coefficients consist of four parts: (i) IC from stellar light, (ii) IC from infrared radiation, 
(iii) IC from CMB, and (iv) prompt emission part. Outside the cylindrical region where the 
GALPROP ISRF is defined, we continue with two components: (i) IC from CMB, and (ii) 
prompt emission. This calculation is numerically the most demanding part of the current 
work. However, it is easy to realize that it can be done only for once and for all: e.g., if 
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one is interested in decaying DM models with a half- life r, then one can simply modify the 
emission coefficients Je by multiplying with a factor of 2mDM/{T{aAv)stdpQ) [98]. Also, it 
is trivial to rescale emission coefficients in order to accommodate various density profiles 
for the Galactic halo. As in Cirelli et al. [83] and in Meade et al. [34], in our calculations 
we also use three distinct forms for the Galactic DM density profile [99, 100, 101]: 



Pq 



exp 



^•0 



1 + 



1~G 



1 + 



!J3 



where the DM density in our local neighborhood pq = p{r = tq 
cm-3^ e.g. [96, 97]. 



Einasto, = 20 kpc, a = 0.17 , 

NFW, rs = 20 kpc , (3.2) 

cored isothermal, = 5 kpc , 

8.5 kpc) = 0.4 GeV 
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Figure 2: The same as Fig. 1 for the decaying DM with half-life r = 5 x 10^^ s. Galactic and 
extragalactic spectra are given with solid and dashed lines, respectively. 

In the following sections, where we compare the models to the Fermi-LAT diffuse 
gamma-ray data, we need to calculate intensities averaged over the annuli defined by ranges 
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of the Galactic b coordinate: 



Aimax 27r 

^Vin<fe<f>.na. ^ L ^ • / d/. / lE{f^,l)dl, (3.3) 

A*min 

where /i = sin(5). In what fohows we will also use an isotropic component of the Galactic 
DM gamma-ray emission. To extract it we find for each photon energy the direction of the 
minimum intensity and use the resulting minimum value as our estimate for the amplitude 
of the isotropic spectrum. As is probably expected, it turns out that the isotropic intensity 
obtained this way practically coincides with an intensity towards b = 0, I = ir, i.e. towards 
the direction directly away from the Galactic center. Even though the spectral shape 
depends on the direction one looks at, the drop in overall amplitude is almost completely 
driven by the rapid falloff of the DM density profile, and so irrespective of the energy, the 
minimum intensity is always practically towards the Galactic anticenter. 

Some example Galactic spectra are shown in Figs. 1 and 2 with solid red lines. In 
each case we have shown the total gamma-ray spectrum along with its decomposition 
into different components. Here the components in the order of increasing energy are: 
(i) IC from CMB, (ii) IC from infrared radiation, (iii) IC from stellar light, (iv) prompt 
gamma-ray emission. In Fig. 1 we have assumed 2/i annihilation channel with a particle 
mass m£)M = 1 TeV and annihilation cross-section {gav) 60 times the standard thermal 
production value of 3 x 10~^^ cm^s~^. The upper and lower left-hand panels show the 
spectra averaged over the Galactic 10 < 6 < 20 and 20 < 6 < 60 regions, respectively. The 
right-hand panels correspond to the 6 > 60 region (top panel) and to the isotropic case 
(bottom panel). In all of the panels the points with errorbars show the measurements of 
the Fermi-L AT [77] . Here all the calculated Galactic models have assumed Einasto density 
profile. Note how the shape of the spectrum changes, particularly notice how the slight 
valley between the IC CMB part and the prompt part of the spectrum gets filled up, as 
one moves from the isotropic case towards the cases that involve more central regions, i.e. 
20 < 6 < 60 and 10 < 6 < 20. The analogous results for the decaying DM with r = 5 x 10^^ 
s half-life are given in Fig. 2. In comparison to the annihilating DM case, where the spectra 
reach E = rriDM, the spectra for the decaying DM extend to as expected. 

3.2 Extragalactic contribution 

Due to the expansion of the space, even in the absence of the absorption, the intensity is 
no longer conserved along the line of sight. In order to calculate the extragalactic contri- 
bution to the gamma-ray intensity due to annihilating/decaying DM we have to include 
this cosmological dimming. Also, no longer can we neglect the absorption of gamma-rays 
once cosmologically large distances are involved. Thus, as in Hiitsi et al. [60], we can write 
for the intensity at redshift z 

oo g 
U.,E) = / cl.- ^,^,^;^^,) (i^) i.(.',£')e>cp [-r(i;,../)] ^ (3.4) 

2 
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Figure 3: The redshift z' where the optical depth for photons reaches unity (i.e., t{E,z,z') — 1 
in Eq. (3.4) ) for several "observer's redshifts": z = 0, 1, 4, 8. Here the energy plotted is the photon 
energy at redshift z. The solid, long-, and short-dashed curves correspond to "no UV", "optimistic 
UV" , and "realistic UV" background cases, respectively. The shaded gray region gives the energy 
range where the Fermi-LAT measurements are done. 



Here the factor {H{z){l + z))~^, where H{z) is the Hubble function, converts the redshift 
interval to proper distance interval, the factor ((l-|-2;)/(l-|-z'))^ accounts for the cosmolog- 
ical dimming of intensities, and the function t{E, z, z') is the optical depth describing the 
absorption between redshifts z and z' . E' = jqi^^-E is the energy of a photon at redshift z', 
assuming it has energy E at redshift z. To calculate the Hubble function we assume a flat 
ACDM "concordance" cosmology with = 0.27 and h = 0.7. In Fig. 3 we have plotted 
redshifts z' where t{E, z, z') = 1 for a wide range of photon energies and for several "ob- 
server's redshifts": z = 0, 1,4,8. Here the energy plotted is the photon energy at redshift 
z, which is the reason why the curves move to right once larger observer's redshift is taken. 
In the lowest energy section of the plot, before the first small plateau, the energy losses are 
dominated by photoionizations. Beyond that, up to the beginning of the large plateau, the 
Compton losses are dominating. At the large flat plateau region the energy loss is domi- 
nated by pair production on matter. In the flnal falling part of the curves the energy losses 
are determined by photon-photon pair production. To calculate the absorption coefficients 
for all of those processes we follow the treatment in Zdziarski & Svensson [102]. The energy 
range relevant for this paper, i.e. the range where the Fermi-LAT measurements are done, 
is given in Fig. 3 by a shaded gray area. Thus, to calculate the optical depth we have to 
take into account the processes of pair production on matter and on ambient photon fields. 
Also, the photon-photon scattering gives some contribution in the region where the flat 
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plateau of Fig. 3 turns over to falling curves. We provide a compact description of how 
we claculate t{E, z, z') in Appendix B. In addition to the CMB photons we have included 
the UV background, which is produced in the lower redshift Universe once the first stars 
start to light up. In Fig. 3 we have shown t{E, z, z') = 1 curves for three distinct models: 
(i) no UV background (solid lines), (ii) "optimistic" UV background (long-dashed lines), 
where for z > 4 we have used the UV background model of Inoue et al. [103], which is 
smoothly joined to lower z estimates of Stecker et al. [104], (iii) "realistic" UV background 
(short-dashed lines), where we have taken into account that recent studies of blazars, e.g. 
[105], suggest significantly lower values for the UV photon densities than estimated in many 
of the previous investigations, and so we reduce the UV background level by an order of 
magnitude for z > 4 and connect this model smoothly to the z = model of Stecker et 
al. [104]. Compared to the recent UV background estimates of Gilmore et al. [106] our 
"realistic" UV model is still somewhat on the "optimistic" side. However, for our study 
these discrepancies are of minor importance, as in reality we have covered a very broad 
range of possibilities: any reasonable model would certainly be somewhere between the 
"no UV" and "optimistic UV" cases. Also, it turns out that for most realistic models the 
extragalactic signal is subdominant, thus leaving the details of the UV background model 
quite irrelevant. 

In comparison to the IC calculations of the previous subsection this time we can sig- 
nificantly speed up our numerics, as a simple Thomson approximation turns out to be fully 
adequate^. Compared to the Galactic case where there are plenty of eV-range photons 
around, here we are fully dominated by very soft target photons of the CMB^. There is 
also a small contribution due to more energetic UV background photons, but those give rise 
to only a negligible contribution to the total IC emissivity. Due to the fact that Thomson 
approximation turns out to be valid, we need to calculate the IC emissivity only for one 
fixed redshift, e.g. for z = 0, and then for other redshifts use simple scaling properties. For 
the case of annihilating DM we can write down for the total emission coefficient 



jr''{z,E) = B{z){l + z) 



1 



l + z 



Jz=0 



E 



l + z 



-prompt / p,N 

+jz=o y.^) 



(3.5) 



where J^Eq ^'^'^ emission coefficients at z = assuming an average DM density p. 

Due to the structure formation the typical value for is not simply but is boosted by a 
large factor B{z), i.e. (p^) = B{z)p^. To calculate the boost function B{z) we use the Halo 
Model [107] of the large-scale structure, which assumes that the global DM density field 
can be modeled as a superposition of DM halos with universal density profiles, and whose 
number density is given by the mass function n(M, z). The details of how we calculate 
B{z) along with several examples can be found in our previous work [60]^*^. Some of the 



See Appendix A for details. 
®Due to the strong cosmological dimming the low redshift gamma-ray intensity is mostly originating 
from a rather local neighborhood only, and thus one does not have to worry about CMB photons becoming 
energetic once z increases. 

^"Note that there is a typo in Eq. (2) of [60]: under the integral sign there should be an additional factor 
of M. 
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results are given at the beginning of Appendix C. The function B{z) is highly sensitive to 
the choice of the concentration- mass relation C{M) for the DM halos. Since the possible 
choices for the C{M) relation already introduce a large uncertainty in the boost function 
B{z) we have chosen to use only one form for the DM density profiles - NFW. B{z) is also 
strongly dependent on the lowest halo mass Afmin- As in our previous study [60] in this 
work we also use two Mmin values: 10~^ and 10~^ Mq, which are quite typical values for 
the WIMP-type DM, e.g. [108, 109]^^. For the concentration-mass relation C(M) we use 
two models: (i) Maccio et al. [95] model, (ii) a simple power-law C{M) model. Thus, in 
total we have four distinct cases for the boost function B{z). 

For the decaying DM, as the energy input is simply proportional to the DM density 
we should replace (1 + z)^ with (1 + z)^ in Eq. (3.5) and also take B{z) = 1, i.e. 



In Fig. 1 we have shown extragalactic spectra for two extreme cases: (i) Myam = IO^^Mq 
along with Maccio et al. C{M) (dashed lines), (ii) Mmin = 10~^Mq and power-law C{M) 
(dotted lines) . At higher energies the curves branch into three distinct possibilities depend- 
ing on the UV model applied: the highest curve corresponds to the "no UV" and the lowest 
to the "optimistic UV" case. The extragalactic spectra for the decaying 2fj, DM model are 
given with dashed lines in Fig. 2. It is worth noting that in comparison to the Galactic 
spectra extragalactic curves have a prominent minimum between the IC and prompt part 
of the spectra, owing to the fact that outside the Galaxy, in the typical intergalactic space, 
we do not have a significant number density of eV-range target photons available, which 
could help in filling up this minimum. 

4. Constraints from Fermi-LAT data 

In this section we present the constraints on leptonically annihilating/decaying DM mod- 
els using the recent gamma-ray data from the Fermi space telescope [77]. We use the 
Fermi-LAT data for the three distinct sky regions defined by the ranges of the Galactic b 
coordinate: 10 < 6 < 20, 20 < 6 < 60, b > 60, along with the estimated isotropic diffuse 
component. In our calculations we use the diffuse signals only, i.e. from the Fermi-LAT 
measurements in the above three sky regions we have subtracted the point source contri- 
bution as given in Abdo et al. [77]. It is clear that the remaining diffuse signal contains 
contributions from many of the known Galactic sources, and thus DM annihilation/decay 
can contribute only some fraction of the total signal. In this paper we have ignored those 
additional complications and have calculated bounds on the models that can be considered 
very conservative: in calculating the constraints we have demanded that in none of the 
energy bins the annihilation/decay signal should exceed the 1-sigma upper values for the 
total diffuse signal as measured by the Fermi-LAT. 

^^As our DM has only leptonic couplings, and A/min is determined by the mass inside the horizon at 
kinematic decoupling, one could arguably reduce Afmin still further, i.e. for the considered scenario our 
values for Mmin are surely on the conservative side. 




(3.6) 
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In a similar manner, the isotropic diffuse component can also have several possible 
"contaminating" Galactic contributions along with extragalactic ones like AGNs^^, struc- 
ture formation shocks [111], emission from starburst galaxies [112] etc. However, it is not 
completely unrealistic that a significant fraction of the isotropic diffuse gamma-ray signal 
might be due to the annihilating/decaying DM. This possibility will be explored in the 
next section, where we show that several DM models are indeed able to provide acceptable 
fits to the observed isotropic diffuse component. 

Our constraints (exclusion zones) on uidm — {cav) plane for five distinct leptonic 
annihilation channels: 2e, 2fi, 2r, 4e, 4//, and for three assumed density profiles for the 
Galactic DM halo: Einasto, NFW, cored isothermal, are given as shaded gray regions in 
Fig. 4. Here the annihilation cross-section {aAv) is given in units of the standard thermal 
production cross-section {cr Av)std — 3 x 10"^^ cm'^s"^. The red ellipses present the 2- and 
1-sigma favored regions of Meade et al. [34] inferred from the measurements of the 
PAMELA, Fermi, and HESS (PFH) experiments. 

In order to determine whether the constraints on Fig. 4 are dominated by Galactic or 
extragalactic component we have plotted in Fig. 5 the ratio of those two for one particular 
case: 2// channel with Einasto density profile. For all the other 14 cases of Fig. 4 the 
picture is qualitatively the same. In general, it turns out that the extragalactic signal can 
dominate only if power-law concentration-mass relation C(M) is assumed. Once arguably 
more realistic Maccio et al. C{M) relation is used the Galactic signal always exceeds the 
extragalactic one. 

In Fig. 4 the strongest bounds are obtained from the isotropic diffuse gamma back- 
ground once the power-law C(M) relation with Mmin = 10~^Mq is assumed. Since in this 
case the extragalactic signal dominates there is a clear sensitivity to the chosen UV back- 
ground model, and of course the strongest bounds correspond to the "no UV" case. The 
10 < 6 < 20 region gives somewhat weaker bounds, and the other regions, i.e. 20 < 6 < 60 
and 5 > 60, which for clarity are not shown in the figure, give constraints between those 
two extremes. Once Maccio et al. C(M) model is used the Galactic signal almost always 
dominates over the extragalactic, and thus there is evidently no sensitivity to the applied 
extragalactic UV model etc. Only if isotropic component is considered the Galactic and 
extragalactic components have a similar magnitude and so there is some sensitivity to the 
choice of the UV model. In contrast to the power-law case the Maccio et al. C{M) gives 
the strongest bounds in the 10 < 6 < 20 region and the weakest ones for the isotropic 
component. Again, the constraints from the 20 < 6 < 60 and 5 > 60 region fall between 
those extremes. 

Keeping in mind that our constraints are very conservative it is clear that in the 
case of the power-law C{M) the PFH favored regions of Meade et al. [34] are rather 
convincingly ruled out. In Fig. 4 we have shown the constraints for the power-law C{M) 
with Mmin = 1O~^M0 only. If Mj^im = 10~^Mq is used instead, the bounds get weaker, 
but the above conclusion stays the same. On the other hand, if Maccio et al. C{M) is 
used, we see that almost always the PFH favored regions survive except for the case of tau 

^■^In a recent paper the AGN contribution is estimated to be around ~ 16% [110]. 
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Figure 4: The constraints on annihilating DM for all five channels: 2e, 2/j,, 2r, 4e, 4/j, and three DM 
density profiles for the Galactic halo: Einasto, NFW, cored isothermal. The gray regions represent 
the exclusion zones for various models as specified in the legend. The red ellipses correspond to 2- 
and 1-sigma favored regions of Meade et al. [34] inferred from the measurements of PAMELA, 
Fermi and HESS experiments. The blue/green ellipses represent the 2- and 1-sigma best-fit regions 
for the isotropic diffuse gamma- ray background for the models without/with Galactic DM halo 
substructure. For the case with substructure we have assumed p°^^ = 15 x Pq. 



lepton final states which are tightly constrained due to the large fraction of the final state 
prompt photons. 
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Figure 5: Ratio of the Galactic to extragalactic signal for the 2/i channel assuming Einasto profile 
for the Galactic DM halo. The extragalactic contribution dominates only if power-law C{M) relation 
is assumed. For the other channels and density profiles the picture is qualitatively very similar. 

The similar constraints for the decaying DM, for the models that are able to fit the 
PFH data are given in Fig. 6. As can be seen from Fig. 5 in the case of decaying DM 
the Galactic signal always exceeds the extragalactic one^'^. Hence the bounds are rather 
insensitive to the choice of the UV model. As there are now less possibilities to consider 
we have shown the constraints for all three Galactic regions, along with the one from the 
isotropic diffuse component. It turns out that the strongest bounds are now obtained from 
the isotropic component and the weakest ones from the 10 < 6 < 20 region. We see that in 
the 2r case the PFH favored region is mostly ruled out, while for the other two cases, 
i.e. 2/i and 4fi, the exclusion zones are still below the PFH ellipses. 



5. Fitting the isotropic diffuse gamma-ray data 

Looking at the model spectra for the isotropic component given on the lower right panels 
of Figs. 1 and 2 one notices that the shapes are quite similar to the observed data, and 
so it is tempting to carry out a fitting exercise. Even though, as mentioned above, there 
are many possible sources like AGNs, emission from starburst galaxies, gamma-rays from 
structure formation shocks etc that might contribute to the isotropic diffuse gamma-ray 
background, it is not completely unrealistic that a large portion of this signal might be due 
to DM. Under that brave assumption we have carried out a fitting procedure using both 
annihilating and decaying DM models. 

^^Note that this result disagrees with the findings of [113], where the authors claim a significantly higher 
contribution from the extragalactic component. 
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Figure 6: The same as Fig. 4 for the decaying DM. Only the models which fit PAMELA, Fermi 
and HESS data are shown. 



As expected, there are indeed regions in the parameter space which provide accept- 
able fits. By looking at Fig. 1 it is clear that in order to obtain an acceptable fit for 
the annihilating DM one has to keep the extragalactic component as small as possible, 
which is achieved by using Maccio et al. C(M) relation. Only the 2t channel, due to 
its large fraction of prompt photons, is not able to provide acceptable fits. For all the 
other channels the best fitting 2- and 1-sigma regions are given as blue ellipses in Figs. 4 
and 6. So far we have assumed that the Galactic DM halo has a smooth density distribu- 
tion characterized by density profiles p(r) as given in Eq. (3.2). However, in case of the 
CDM cosmologies we know that structure formation proceeds according to the hierarchical 
bottom-up scenario, and thus one would expect significant amount of substructure (subha- 
los) inside our Galaxy's DM halo. As discussed e.g. in Martinez et al. [109] , depending on 
a particular substructure model, these subhalos might potentially boost the smooth halo 
annihilation signal by a factor of ~ 1 — 20. In our calculations we will assume a simple 
model for the halo substructure, where the subhalos follow the main halo's density profile 
and also that subhalo mass distribution function along with concentration-mass relation is 
independent of the location inside the main halo. Under those assumptions, as shown in 
Appendix C, the gamma-ray emissivity due to halo substructure scales as oc /^(r), i.e. sim- 
ilar to the case with decaying DM. As shown in Appendix C under the above simplifying 
assumptions the specific details of the substructure model can be absorbed into a single 
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Figure 7: Best fit annihilating DM models (with Einasto density profile) for the isotropic diffuse 
gamma-ray background (upper left panel). The solid and dashed lines represent the models without 
and with the Galactic substructure, respectively. For the case with substructure we have assumed 
Psuh — 15 X /90. The averaged spectra of the same models in several Galactic regions are also shown: 
6 > 60 (upper right panel), 20 < 6 < 60 (lower left panel), 10 < 6 < 20 (lower right panel). Dotted 
lines show the standard Galactic diffuse background (IC + bremsstrahlung + Tr^-decay) model 
44_500180 from Strong et al. [114]. NOTE: The fit is performed only for the isotropic diffuse 
component (upper left panel). For the other regions the minimal requirement for acceptability is 
that models should not exceed the data points. 

location-independent parameter pl^^, allowing the calculations to proceed practically the 
same way as for the case without substructure, only p^(r) needs to be replaced by p'^^-^^- p{r). 
If one takes p^^^^ = P& = 0.4 GeV cm^^ then the signal from substructures turns out to 
be approximately equal to the main halo contribution, and in loose terms one can say 
that the substructure boost (1 + substructure signal / main halo signal) is approximately 
equal to 2. In order to lower the blue best-fit regions of Fig. 4 close to the PFH favored 
red ellipses one has to have substructure boost factors of order ~ 10. In Fig. 4 the green 
regions correspond to the case where p^^-^ = 15 x pQ. 
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The best fitting anniliilating DM model spectra for tlie 2e, 2/i, 4e, and cliannels are 
given on tlie upper left-hand panel of Fig. 7. The solid and dashed lines represent the models 
without and with the Galactic substructure, respectively. For the case with substructure 
we have assumed p'^^-^ = 15 x pQ. In the other panels we have plotted the spectra for the 
same models in three Galactic regions: 6 > 60 (upper right panel), 20 < 6 < 60 (lower 
left panel), and 10 < 6 < 20 (lower right panel). We see that although the models can 
nicely fit the isotropic diffuse component, the ones without substructure clearly exceed the 
observational data once more central regions of Galaxy are included, i.e. 20 < 6 < 60 
and 10 < 6 < 20 regions. Once substructure is included those problems with central 
regions start to disappear. In reality, the substructure does not need to simply follow 
the density profile of the main halo, but in the central regions the destruction probability 
for subhalos might significantly increase. This helps in lowering the central signal with 
respect to the contribution from the outskirts of the halo, and thus making the gamma-ray 
emission profile shallower. Note that the profile can also be made shallower if one allows 
the intermediate particle (j) (see Section 2) to have a sufficiently long lifetime, so that it 
can move out of the initial place of production before decaying, which would effectively 
correspond to the smoothing of the gamma-ray emission profile^^. Following this line of 
thought it is not hard to imagine the case where the problems with the central regions 
in Fig. 7 disappear altogether. Thus, it seems quite possible to have annihilating DM 
models, which: (i) are compatible with the PFH data, and (ii) provide acceptable fits 
to the Fermi-LAT isotropic diffuse background, without being in clear conflict with the 
other Fermi-LAT gamma-ray measurements. As the modern N-body simulations are only 
capable of resolving DM substructures that are several orders of magnitude larger than 
the mass scales relevant for the current problem, one is necessarily limited to this type of 
qualitative arguments only. 

The analogous results for the decaying DM are given in Fig. 8. Here in contrast to the 
annihilation case (without substructure) we see that the best fitting models of the upper 
left-hand panel do not exceed the observational data points in other regions. However, the 
blue best-fit regions in Fig. 6 are quite far from the PFH ellipses. 

6. Inclusion of the standard Galactic background 

So far we have considered only DM-induced gamma-ray signal. The minimum requirement 
for the DM models to be valid is not to exceed the observational data points in Figs. 7 
and 8. However, it is clear that diffuse gamma-ray background has other contributions. 
Here we consider the standard Galactic background model 44_500180 from Strong et al. 
[114]. In Figs. 7 and 8 with dotted lines we have plotted this standard Galactic diffuse 
background model, which consists of IC, bremsstrahlung, and vr'^-decay components. In 
Figs. 9 and 10 we have added this background to the DM-induced signal. 

^*0f course, the acceptable range of lifetimes is limited, as in the most radical case, where the lifetime of 
<j> gets very large, the Galactic signal becomes very diluted and the typical IC spectrum would be generated 
on extragalactic photon fields, instead. 
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Figure 8: The same as Fig. 7 for decaying DM. 

It is quite remarkable that the models, which provide the best fits to the isotropic 
diffuse data, automatically get quite close to the observed data points in other Galactic 
regions, once the standard background model of Strong et al. [114] is included. It is 
evident that by a relatively modest tuning of the Galactic background model satisfactory 
fits in all of the regions should be possible. The only exception is the annihilating DM 
case for the region 10 < 6 < 20, where at the highest energies the models are above the 
observational data points, even before the standard Galactic background is added. But 
this, as explained earlier, could in principle be alleviated by, e.g., increasing the dominance 
of the substructure at the outer halo regions, this way effectively reducing the steepness of 
the resulting gamma-ray profile. 

7. Summary and discussion 

In this paper we have used the gamma-ray data from the Fermi-LAT space telescope to 
place constraints on leptonically annihilating and decaying DM models. In particular, we 
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Figure 9: The same as Fig. 7 with the standard Galactic background added to the DM signaL 



have used the diffuse gamma-ray data for the Galactic regions 6 > 60, 20 < 6 < 60, and 
10 < 6 < 20, along with the isotropic diffuse component as estimated in Abdo et al. [77]. 
Thus, in our calculations we do not consider the very central regions of the Galaxy, as the 
closest we get is 10 degrees from the center. This new Fermi data has been also used in 
some of the other analogous studies. In contrast to the present study, where we calculate 
both prompt and IC parts of the gamma-ray spectra, the authors of [89] have focused only 
on the prompt emission part. In Abdo et al. [90] the authors have paid a significant amount 
of attention to the extragalactic contribution, but their Galactic signal does not seem to be 
very realistic, as the IC and prompt components of the spectra are clearly separated, which 
hints that their model for the ISRF used in the IC calculations, has probably neglected 
the stellar radiation and the infrared radiation from the Galactic dust. In comparison 
to our paper, were we have investigated both decaying and annihilating DM, the above 
studies consider only the annihilating case. The preliminary Fermi data has been used to 
constrain DM models in Meade et al. [34], in Cirelli et al. [83], and in Papucci & Strumia 



-19- 



0.01 



CO 

0.001 

E 
o 

> 

CD 



in 




0.0001 
0.01 



LU 
LU 



0.0001 



decay, Einasto profile 



decay, Einasto profile 



Fermi isotropic diffuse i 
2fi channel, rripf^ = 1 480 GeV, t = 5.99 x 1 0^^s 
4(1 channel, m^f^ = 1 230 GeV, i = 2.43 x 1 s 



isotropic 



20 < b < 60 




Fermi 20 < b < 60 region i 




b>60 



1 < b < 20 



100 1000 10000 

E [MeV] 



100000 



Fermi b > 60 region i 



Fermi 10 < b < 20 region i 




1000 10000 

E [MeV] 



100000 



Figure 10: The same as Fig. 8 with the standard Galactic background added to the DM signal. 

[82]. Compared to our previous study [60] where we focused only on extragalactic signals, 
this time we have included the contribution from our own Galactic DM halo. The most 
important difference between our analyses and the previous papers is that we have also 
performed fits of the diffuse gamma-ray data and compared the results with the preferred 
DM parameters explaining the PFH anomalies. 
Here are our main conclusions: 

• For the annihilating DM models the extragalactic signal can dominate only if power- 
law concentration-mass relation C(M) is assumed. Once arguably more realistic 
Maccio et al. C{M) relation is used the Galactic signal always exceeds the extra- 
galactic one. Due to those reasons our calculated bounds on annihilating DM have 
significant dependence on the assumed extragalactic UV background model (which 
strongly effects the level of prompt emission) only in the case of the power-law C{M) 
relation. In the case of decaying DM, where the signal scales proportionally to the 
DM density, and thus is not very sensitive to the details of the structure formation 
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model, the Galactic signal is always stronger than the extragalactic one. As a re- 
sult the calculated bounds in this case are rather insensitive to the choice of the 
extragalactic UV model. 

• The derived constraints are very conservative, as the only thing we have required 
is that none of the models should exceed the measured diffuse signals in any of the 
observed energy bins. It is clear that in reality those diffuse signals can have various 
Galactic and extragalactic contributions other than the potential contribution from 
the annihilating/decaying DM. Our constraints for the annihilating and decaying DM 
models are given in Figs. 4 and 6, respectively. The red elliptical regions in those 
figures show the PFH favored regions of Meade et al. [34] . It turns out that for 
the power- law C{M) relation the PFH favored regions are quite convincingly ruled 
out. On the other hand, if Maccio et al. C{M) is used, we find that almost always 
the PFH favored regions easily survive except for the tau final states that produce too 
many prompt gamma-rays. For the decaying DM only the 2r channel is significantly 
constrained, while for the other two cases, i.e. 2/i and 4/i, the exclusion zones are 
still below the PFH ellipses. 

• Once realistic models for the ISRF are used with all the relevant components: (i) 
stellar light, (ii) infrared radiation from dust, (iii) CMB, it turns out that the re- 
sulting gamma-ray spectra from the Galactic annihilating/decaying DM halo can be 
remarkably close to a seemingly featureless power law over the relevant energy range, 
as seen in Figs. 1 and 2. Thus, there might be a possibility that those models provide 
an acceptable fit to the isotropic diffuse gamma-ray data if extragalactic contribution 
can be kept low, which is surely true for the Maccio et al. C{M) relation. Indeed, this 
turns out to be the case: in Figs. 4 and 6 the blue ellipses show the 2- and 1-sigma 
best-fit regions for the scenario neglecting Galactic DM halo substructure. Only the 
r lepton final states fail to give reasonable fits to the data due to hard prompt pho- 
tons. However, there is no overlap between the gamma-ray and PFH best-fit regions. 
Several best fitting spectra corresponding to this case are plotted with solid lines on 
the upper left-hand panels of Figs. 7 and 8 for the annihilating and decaying DM, 
respectively. It is clear from Fig. 7 that the best fitting annihilating models are in 
conflict with other measurements from other regions, especially with the ones from 
the 10 < 6 < 20 Galactic region. However, in the case of the decaying DM there is 
no immediate conflict (i.e. the model spectra do not exceed the data points) with 
the other measurements, as seen from Fig. 8. 

• All the above results did not include the possible contribution from the Galactic 
DM substructure. We have devised a simple model for the halo substructure (see 
Appendix C), where all the specific model details can be absorbed into a single 
effective substructure density parameter /OgJ^^. The inclusion of substructure has 
a twofold effect: (i) it leads to the increase of the total amplitude of gamma-ray 
emissivity, (ii) it makes the spatial emission profiles shallower compared to the smooth 
halo case. Due to those reasons the inclusion of substructure lowers the blue best- 
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fit regions of Fig. 4, and at tlie same time helps to reduce the conflict with the 
measurements of other Galactic regions seen on the lower panels of Fig. 7. The 
green best-fit regions of Fig. 4 and the dashed lines of Fig. 7 correspond to the 
case with p'^^^ = pq = 0.4 GeV cm~^, which in loose terms represents the case 
with substructure boost factor ~ 15. Those levels of substructure boosts are indeed 
possible, as discussed, e.g., in Martinez et al. [109]. If we give up our simplifying 
assumption that substructure follows the density profile of the main halo, and allow 
an increased subhalo destruction rate in the center, it should be possible to make the 
spatial emission profiles even shallower, and thus potentially remove the remaining 
small conflict on the lower right-hand panel of Fig. 7. The emission profiles can also 
be made shallower if one allows the intermediate particle (j) (see Section 2) to have a 
sufficiently long lifetime, so that it can move significantly away from the initial place 
of birth before decaying, which would effectively reduce the steepness of the resulting 
gamma-ray profile. 

It is clear that diffuse gamma-ray background has other contributions. In this work we 
consider the standard Galactic background model 44_500180 from Strong et al. [114]. 
In Figs. 9 and 10 we have added this background to the DM-induced signal. It is quite 
remarkable that the models, which provide the best fits to the isotropic diffuse data, 
automatically get quite close to the observed data points in other Galactic regions, 
once the standard background model of Strong et al. [114] is included. It is evident 
that by a relatively modest tuning of the Galactic background model satisfactory fits 
in all of the regions should be possible. 

We conclude that it seems quite realistic to have annihilating DM models, which 
simultaneously: (i) are compatible with the PFH data, and (ii) provide acceptable 
fits to the Fermi-LAT isotropic diffuse background, without being in clear confiict 
with the other Fermi-LAT gamma-ray measurements. The preferred DM annihilation 
final states consistent with this analysis and the PFH fits are 2p, 4/i and 4e. 

• In our Galactic IC calculations we have neglected the effect of electron diffusion, i.e. in 
our treatment the IC spectrum is produced "on spot" . For more precise treatment one 
should solve the diffusion-loss equation for electrons to determine the electron energy 
distribution, which is used as an input in IC calculations. Some estimates of Papucci 
& Strumia [82] have shown that by including this diffusion effect the IC signal gets 
reduced by a factor of around 2—3 depending on the Galactic region. However, the full 
calculation of electron diffusion is quite complicated with rather large uncertainties 
due to, e.g., the knowledge of the magnetic fields. Even in our simple treatment there 
is a large uncertainty due to the assumed ISRF model. We have used the latest ISRF 
model from the GALPROP code, which is based on results from Porter &: Strong 
[86] . However, in their study the authors do not provide any errors for the estimated 
ISRF. It is clear that their results are highly model dependent, as from the direct 
measurements we can only get the ISRF in our location, while the values at other 
Galactic locations can only be obtained by solving the radiative transfer problem, 
which clearly involves several model-dependent assumptions. According to Strong 
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et al. [114] a factor of two uncertainty for the ISRF is quite possible. Due to those 
reasons we feel that the accurate treatment for the diffusion problem can probably 
wait before there are better models with realistic uncertainties available for the ISRF. 
Unfortunately, except for rare examples of [86] and [115], the problem of determining 
the ISRF as a function of Galactic position has not gained much attention so far. 

• In Introduction we referred to the other constraints on annihilating or decaying DM. 
Compared to the CMB constraints (see e.g. [59, 60, 61]) the bounds derived in this 
paper can be stronger or weaker, depending on what type of concentration-mass rela- 
tion C(M) one assumes. If Maccio et al. [95] C(M) is used the bounds are generally 
weaker, while the power- law C(M) leads to stronger constraints than obtainable from 
the CMB measurements. The CMB constraints, being mostly determined by the free 
electron fraction at redshifts 100 ^ z < 1000 [60], are quite insensitive to the un- 
certainties of the nonlinear structure formation and other complicated low-redshift 
astrophysical phenomena. In that respect they can be considered significantly more 
robust than the diffuse gamma-ray bounds derived in this paper. 

The constraints from the Galactic center radio and gamma-ray measurements are 
generally stronger than ours, however those are hampered by the obscurity of the 
DM density profile, ISRF and magnetic field in the very central region (see, e.g. 
[28, 29, 30, 47, 48]). For example, there is no clear understanding on how much various 
physical processes involving baryonic component can impact on the the central cusp 
of the DM profile (see e.g. [116] for a recent review). 

The constraints derived from the Galactic dwarf spheroidal galaxies and the nearby 
galaxy clusters are comparable or slightly weaker than our restrictions [42, 41, 46]. 

Note added 

When the research presented in this work was completed, Ref. [117] appeared that also 
performs fits of the cosmic ray and gamma-ray data. Although their approach and aim 
of the paper is different form ours, the main conclusions that only DM annihilations can 
potentially provide a satisfactory fit agree with each other. 

Acknowledgments 

We thank Alessandro Strumia, Marco Cirelli and Mario Kadastik for discussions and our 
Referee for useful comments and suggestions. This work was supported by the ESF Grants 
8499, 8090, 8005, 7146, by SF0690030s09 and by EU FP7-INFRA-2007-1.2.3 contract No 
223807. GH acknowledges the support from the EstSpacE fellowship. 

A. Emission coefficients for DM annihilation and decay 

The total gamma-ray emissivity j^; [eV cm~^s^"'^eV^^] for the leptonically annihilating/decaying 
DM models considered in this paper can be given as a sum of IC and prompt parts 

jE{E)=jf{E)+f^""'^\E), (A.l) 
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where the prompt emission part is trivial to calculate once the input photon spectra from 
PYTHIA simulations are available. In this appendix we focus on the IC calculations. In 
what follows we assume that the IC spectrum is generated "on spot", allowing us to use 
one fixed value for the DM density, p, as well as for the number density of target photons 
with energy E . . .E + dE, nE{E). The IC emissivity can then be expressed as 



jf{E) 



{<^A'v) p2 £qj. annihilating DM 



p/t 



for decaying DM 



dE, 
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'dE, 



(A.2) 



where which is calculated with PYTHIA, is the spectrum of produced per annihi- 
lation/decay. The functions Ki{E,Ee) and K2{E(,) are given as: 
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B. Absorption coefficients and optical depth for gamma photons 

As described in the main text, the processes relevant for the absorption of gamma-rays in 
the energy range probed by the Fermi-LAT are: (i) pair production on matter, (ii) photon- 
photon scattering, and (iii) pair production on ambient photon fields. At those energies 
by far the dominating one amongst the three is the pair production on photon fields. Here 
we provide a brief summary of results needed to calculate the optical depth t{E, z, z') in 
Eq. (3.4) for gamma-ray photons between redshifts z and z' . For for further details see 
[102] and [118]. In what follows we express the gamma photon energy in units of electron 
rest mass, i.e. e = — 
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The absorption coefficient at redsfiift z for pair production on neutral matter witli 
mass fractions of hydrogen X = 0.75 and helium Y = 0.25 can be approximated as 

amat-pair(e,2) ^ ao(l + ^)^ln (^^^^^^ , e>6, (B.l) 
«o ^ S.Sn^ayrg 2.05 x lO'^ {J^ {J^) ^pc"! , (B.2) 
where the average electron number density at z = is 



2 



n° 2.17 X 10-^ (-%-)( ) ■ (B-S) 



0.045 J V0.7 

Here Qfe is the density parameter for baryons, h is the Hubble parameter, ctf = 
7.29735 X 10~^ is the fine structure constant, and vq = 2.8179 x 10"^^ cm the classical 
electron radius. 

However, below redshift z ~ 6 the Universe is known to be reionized [119]. The ab- 
sorption coefficient for pair production on fully ionized matter can be approximately 
given as 



aionmat-pair(e, z) ~ ao{l + Z 



,3 



■ , , 109 
H2e) - - 



e > 1 , (B.4) 



ao - f n^.ro^ ^ 2.58 x 10^^ (^) (A) " Mpc^ . (B.5) 



^ ~ 3.23 X 10^-^1 ° Mpc"i , (B.7) 

455625 Ac V 2-725 ^ ' ^ ^ 



The absorption coefficient for photon-photon scattering at redshift z can be given by 
a77-scat(e, z) = ao(l + zfe^ , (B.6) 

ao -- 

where 0o = is the CMB temperature at z = in electron rest mass units, and 
Ac = 2.4263 X 10^^*^ cm is the electron Compton wavelength. 

For the photon-photon pair production, in addition to the soft CMB photons, we 
have to include photons of the extragalactic UV background, which are produced 
in the low redshift Universe once the first stars start to light up^^. As there are no 
simple analytical forms available for the UV background spectrum, we have to use full 
numerical treatment for the photon-photon pair production process. The absorption 
coefficient at redshift z for photon-photon pair production can be given as 



An J ^'''^''^^(eg)2' 



a77-pair(e,^;) = / den{e,z)j-^ , (B. 



^^The influence of those UV photons on the previously discussed photon-photon scattering process, which 
itself is a weak effect, is completely negligible. 
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(l)(v) = — Inw -ln^w + 21n^(l+w) + 
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+ 4Li2 -— --, (B.9) 
\l + w J 3 
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where n(e, z) = ncMB(e; -z) + '^uv(e) is the number density of target photons, and 
Li2 is the dilogarithm function. 

Now we can write for the total absorption coefficient 

/ N J "mat-pair (e,^;) if 6 < Z < 1000 1 f \ , f \ 

[ aionmat--pair(e, z) \{ Z < Q J 

(B.ll) 

and for the optical depth t{E, z, z') between redshifts z and z' 

z' 

r{E, z, z') = c [ dz 'I , (B.12) 

z 

where e = j-}^ ^ o . To calculate the Hubble function H(z) we assume a flat ACDM 

l + Z rUeC^ ^ ' 

model with = 0.27 and h = 0.7. 



C. Halo substructure: a simple model 

In our previous paper [60] we showed that within the Halo Model of the large-scale structure 
the typical value of at redshift z gets boosted by a factor B{z) above the average DM 
density squared, i.e. {p'^{z)) = B{z)fP{z), where 

oo 

B{z) = l + -^ J dMM—{M,z)fc[C{M,z)]. (C.l) 

Here Ac is the overdensity (with respect to the mean DM density) used for defining halo 
masses in simulations, with a typical value of Ac — 200. Mmm is the lower cutoff for the 
halo mass, ^ is the halo mass function, and the function fc, which depends on the halo 
concentration-mass relation, is given for the NFW halos as 

(C.2) 

This calculation looked at the contribution from the ensemble of halos, but did not 
include the substructure inside each of these. It is clear that similar calculation can be 
repeated for one particular halo plus an ensemble of subhalos inside it. One just has to 
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replace the halo mass function with subhalo mass function Thus, for the substruc- 

ture boost factor as a function of distance from the Galactic center we can write (we do not 
include +1 as this corresponds to the contribution from the mean smooth DM component) 

oo 

AT dn 
Bsnhir) = ^ / dM,^^M,^k— (M3ub,r)/c[C(Msub,r)] . (C.3) 

p J G.IV1 sub 

-^-^sub,min 

Now we make assumptions that the subhalo mass distribution function /sub(-^sub) 
(J dMsub /sub (-^sub) = 1) does not depend on Galactocentric radius r, and that the spatial 
distribution of subclumps follows the density profile p(r) of the main halo. We also as- 
sume that the concentration-mass relation C(Msub,'^) does not depend on r. Under those 
simplifying assumptions we can express the subhalo mass function via 

dn , s I'm fsuhiMsnh) f s a\ 
{Msuh,r) = — p{r) . (C.4) 



dMsub' ' Msub 

Here /m is the ratio of the total substructure mass and the main halo mass. For the 
substructure boost factor we can then write 

oo 

Bsubir) = /a/A,^ / dMsub/sub(Msub)/c [C(Msub)] • (C.5) 
P J 



sub, mill 



Here the only r dependent term is p{r). 

In our calculations we need (Psub('^))' which can be be given by 

ipl^ir)) = B,^y,{r)p' = • Pir) , (C.6) 

where the effective substructure density 

oo 

pfb = ImA^p^ I dM3ub/sub(Msub)/c [C(M,ub)] • (C.7) 

-^^sub,min 

Thus, under the above simplifying assumptions all the model details can be absorbed 
into one r independent parameter p^yb- If we take p^^-^ = P& = 0.4 GeV cm~^ then the 
gamma-ray signal from the main halo and from substructures turn out to be of very similar 
magnitude. However, note that the spatial profiles of the signals are quite different: the 
gamma-ray emissivity profile from substructures is similar to the case of the decaying DM, 
as in both cases j£;(r) oc p{r). 

We conclude that the inclusion of DM substructure will generally: (i) increase the 
total amplitude of gamma-ray emissivity, (ii) make the spatial emission profiles shallower 
compared to the smooth halo case. 
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